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Measurement of D^ decays to K^tt^ and K'^n^ at Belle 



The Belle Collaboration 
Abstract 

We present a preliminary measurement of the ratio of D^ decay rates into K^tt^ and K^-k^ 
final states. This ratio can be used to disentangle the Cabibbo favored D^ -^ K tt^ and doubly 
Cabibbo suppressed D^ — > K^tt^ amplitudes, and contributes to the important goal of constraining 
the strong phase 6k-k between D^ -^ K~t:'^ and D^ -^ K'^Tr~. The measurement is based on a 
large data set accumulated by the Belle detector at the KEKB e+e~ collider, with K^ candidates 
reconstructed using hadronic clusters in the KLM (Er£ and fi identification system), together with 
a D^ mass constraint and D* tag. 
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I. INTRODUCTION 

The search for D^ — D mixing has long been one of the most interesting projects in charm 
physics. Recent pubUcations by the FOCUS [jj and CLEO |j] collaborations have stimulated a 
series of new measurements, including an updated ycp measurement by the Belle collaboration [3], 
reported at the summer conferences. However, the experimental situation is still not clear. 

An important quantity in the interpretation of D^ — D mixing results is the phase difference 
^Kn between the Cabibbo favored (CF) and doubly Cabibbo suppressed (DCS) neutral D meson 
decays. As has been shown in |j|, new information on Skw may be obtained by measuring the 
asymmetry between the decay rates of D^ into Kgir^ and -ff^vr", where the effect may be as large 
as 0{tan^ 9c), i.e. at the 5% level. 

In this paper we present the first measurement of the D^ -^ Kf^ir^/D^ -^ K^tt^ decay rate 
asymmetry, using the Belle detector at the KEKB e"'"e~ collider. In order to cancel out most of 
the systematic effects of the detector, we extract the K^/Kg relative detection efficiency from the 
ratio of K^7r~^Tr~ and i^^7r"'"7r~ modes via the decay D^ — > K*~7r~^. The presence of K*~ -^ K^'k~ 
in the decay chain ensures equal rates of K^ and K^ in this case. 

Throughout this paper, the inclusion of CP-conjugate processes is implied. In particular, we 
combine D^ -^ K^^tt^ with D — > K^n^, and so on. We thus treat the X^ and Kg as CP-eigenstates 
(the error is negligible) and assume that direct CP violation in Z) — > Kn and Kit decays, if present, 
may be ignored at our present sensitivity. 

II. Kl RECONSTRUCTION TECHNIQUE 

The key point in the analysis of final states containing if £ is the reconstruction of its momentum. 
The technique presented below makes use of the KLM {Kj^ and /i identification) subsystem of the 
Belle detector (5|, which consists of several layers of iron absorber interspaced with resistive plate 
chambers. The system enables one to reconstruct muon tracks as well as hadronic showers. Those 
showers which do not have a matching charged track are assumed to be produced by if£, with 
the position of the shower giving information about the flight direction of the K^. The direction 
resolution is improved if the K^ deposits energy in the electromagnetic calorimeter in front of 
the KLM, in which case the position of the calorimeter cluster is used to determine the K^ flight 
direction. 

To extract the magnitude of the K^ momentum one needs to employ additional constraints. 
In particular to reconstruct D^ decays into Kf^X where X is a fully reconstructed system, one 
can apply a D^ mass constraint and solve the resulting 4-momentum equation with respect to the 
magnitude of the Kj^ momentum, and then exploit the D* — > D^tt decay to tag the signal. 

To monitor the performance of the method we utilize the same D^ decay modes but with Kg 
instead of if£ in the final state. We use the Kg flight direction reconstructed by its decay into 
7r"'"7r~, and then apply the same procedure as for if£: henceforth such Kg are referred to as -ff£ 
(pseudo K^). 

Since the equation for the magnitude of the K^ momentum is quadratic there may be up to 
two solutions. However one of those, even if present, is non-physical in more than 90% of events. 
For simplicity, in this analysis we use only the one with the larger momentum. 



III. DATA SET AND EVENT SELECTION 

The data sample used for this analysis amounts in total to an integrated luminosity of 23.6 fb~ . 

We select D* candidates with reconstructed scaled momentum Xp > 0.6, which rejects D* 
from i?-meson decays and supresses the combinatoric background. In addition to the standard 
reconstruction procedure for tt^, tt*^, and Kg, the latter are required to satisfy the following quality 
criteria: 

• the K*^ —5- 7:~^ir~ vertex should be separated from the interaction point in the plane perpen- 
dicular to the beam axis by more than 500 /xm; 

• the distance along the beam axis between the tt tracks at the Kg vertex must be less than 
1 cm; 

• the cosine of the angle between the assumed Kg flight path from the interaction point to the 
decay vertex and the reconstructed Kg momentum in the transverse projection with respect 
to the beam must be more than 0.95. 

When reconstructing Kf^, there is a potential background due to KLM clusters caused by unrecon- 
structed charged particles, in particular when these particles do not originate from the interaction 
region (e.g. K and vr decays in flight). To reject such cases we veto any KLM cluster matched 
to a calorimeter cluster with an energy between 0.15 and 0.3 GeV, corresponding to minimum 
ionization. 

The dominant part of the background is due to the pickup of arbitrary soft pions when forming a 
D^ candidate. A characteristic feature of such combinations is that the D^ mass is accommodated 
by assigning a large momentum to the K^, which then accounts for most of the D^ momentum. 
Therefore the K^ flight angle with respect to the D^ boost (hereafter 6dk) for this background 
tends to peak in the forward direction (Fig. ^) while for the signal the underlying distribution is 
isotropic, with a slight tilt to the forward direction in the measured distribution (Fig. ^) due to 
detector efficiency. 
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FIG. 1: A'° flight angle with respect to the D^ boost for D^ -^ K^tt^ 



The sharp peak around cosOdk = — 1 is due to the degenerate case when a random vr*^ or 
7r~^7r" system can form a "good" D^ with a K^ (almost) at rest in the lab: to exclude these cases. 



we require cos O^k > —0.95. The upper cut on this quantity needs optimization which has been 
performed using signal MC and data in the M(D*) sideband: we require cosOdk < 0.2 for ah four 
modes. 

The invariant mass of the K^'k~ combination in the K^t:~^'k~ mode is required to he within 
50 MeV/c^ of the nominal Er*(892) mass. In addition the invariant mass of the Tr"*"?:" pair is 
required to be less than 0.7 GeV/c^ to make the signal and the control modes kinematically 
similar, and to suppress the contribution from K^p decays. 

Since the kinematics in the K^ and the corresponding K^ modes is the same we can use the 
fully controlled K^-k^ and K^tttt modes to select a region in the K^ lab momentum where the 
efficiency of the kinematical cuts, which is independent of the K^ flavor, is the same (see Fig. EJ. 
Hence we can average over this region and all the efficiencies including that of K^ reconstruction 
cancel out. We select the K^ lab momentum range from 0.6 GeV/c to 2.5 GeV/c for this analysis. 

IV. DETERMINATION OF Kl AND K^ RATE RATIO 

The resulting D* mass plots for the four D^ decay modes under study are shown in Figure 01 
For the modes including K^ these are mass difference plots in the traditional reconstruction with a 
very loose cut on the D^ mass (within 100 MeV from the nominal value). The distributions for the 
K^ modes have been fit to a sum of two gaussians representing the signal and its tails and a first 
order polynomial multiplied by a square root threshold factor to describe the background, while 
for the signal in the K^ modes a single gaussian was used. The central value of all the gaussians 
was fixed at the nominal charged D* mass; all other parameters were allowed to float. 

The yields obtained from the fit are 

N{K^sT^^) = 4715 ±91, 

N{k1tt°) = 1839 ±101, 

iV((i^^7r)7r) = 2524 ±77, 

N{{Kl7r)-K) = 1119 ±84. 

They are used to calculate the ratio of branching fractions of D^ -^ K^ir^ and D^ -^ K'^tt^ 

B{D^ -. i^£vr°) _ N{Kln^)/N{Kln^) _ 

B{D^ ^ i^o^O) " iV((K0vr)7r)/iV((if 07r)7r) " ^'^^ U.uy^statj 

V. EVALUATION OF SYSTEMATIC ERRORS 

Due to the similarities in the behavior of the signal K^-k^ and control (i^'^7r)7r modes most of 
the systematic errors cancel out. However some steps in the analysis are different and thus may 
have different systematic effect on the result. 

We consider the following potential sources of systematic errors which may not cancel: 

• the effect of the residual difference in K^ lab momentum spectra, which can lead to a bias 
in the result if the K^ reconstruction efficiency strongly depends on momentum; 

• systematic errors due to the imperfect fitting function. 
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FIG. 2: K lab momentum spectra for (a) Kgii and (b) (i^^vrjvr; (c) shows the ratio of (a) to (b) 
as a function of momentum. 



In order to estimate the bias induced by the residual difference in K^ lab momentum spectra, 
an artificial momentum-dependent "efficiency" was introduced in K^t:' and (i('^7r)7r modes by 
applying a weight depending on the K^ lab momentum. The weight changed linearly from at 
0.6 GeV/c to 1 at 2.5 GeV/c, representing the worst possible case. The resulting change in the 
ratio of yields was less than 3%. 

Systematic errors arising from imperfections of the fitting model were evaluated by varying the 
functions describing the signal and the background: in particular we employed a double gaussian 
to describe the signal in K^ modes as well as in K^, and a second order polynomial multiplied by 
the square root threshold factor to fit the background. The yield ratio remained stable within 6% 
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FIG. 3: D* mass plots for the specified D^ decay modes 

under all variations of the fitting model that were studied. 

In addition, the stability of the fit against variation of the background level was studied by 
simultaneously varying the cut on cos 6dk in all the four modes: we estimate the error on the yield 
ratio due to this source to be 8%. 

Adding these sources of error in quadrature, we estimate the total systematic error to be 10%. 
Since this error is dominated by the difficulty of parametrizing the background, there is potential 
for improvement if the shape of the background can be simplified. One technique currently being 
studied is the replacement of the M{D^) constraint with a constraint on M'^{D*) — M'^{D^) when 
forming the K^ momentum: the signal yield is then extracted from the resulting M{D^) distribu- 
tion. As there is no natural threshold in this case, we expect the yield to be relatively stable under 
variations of cuts, fit functions and so on. 

VI. CONCLUSION 



In summary, we have demonstrated the viability of this method of measuring the ratio of D^ 
decay rates into K^-k^ and K^tt^ . Our preliminary measurement of this ratio is 



B{DO 



Khf) 



0.88 ± 0.09(stat) ± 0.09(syst); 



expressing the result in terms of tlie rate asymmetry defined in [J] , we find 

A ^ r(i^o _ 4vr0) + rpo ^ 4^0) = 0-06 ± 0.05(stat) ± 0.05(syst), 

which is consistent with unity. At the current level of precision we are therefore not able to place 
any strong constraint on the parameters of D ^ Ktt decays, and so the strong phase 6ktt- However 
the statistical error will soon improve as more data is accumulated by the Belle detector; changes 
of technique to reduce the systematic error are also being actively studied. 
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